Mon. Not. R. A.stron. Soc. 000,[T]l6](2010) Printed 4 June 2010 (MN L?TbX style file v2.2) 



High- velocity feature of the class I methanol maser in G309.38— 0.13 
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ABSTRACT 

The Australia Telescope Compact Array (ATCA) has been used to map class I methanol 
masers at 36 and 44 GHz in G309.38— 0.13. Maser spots are found at nine locations in an 
area of 50"x30", with both transitions reliably detected at only two locations. The brightest 
spot is associated with shocked gas traced by 4.5-/im emission. The data allowed us to make 
a serendipitous discovery of a high-velocity 36-GHz spectral feature, which is blue-shifted by 
about 30 km s^^ from the peak velocity at this frequency, but spatially located close to (within 
a few arcseconds of) the brightest maser spot. We interpret this as indicating an outflow par- 
allel to the line of sight. Such a high velocity spread of maser features, which has not been 
previously reported in the class I methanol masers associated with a single molecular cloud, 
suggests that the outflow most likely interacts with a moving parcel of gas. 
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1 INTRODUCTION 

Methanol masers are commonly found in massive star-forming re- 
gions, with more than twenty different centimetre a nd millimetre 
wave length masing transitions discovered to date (e.g.. lMuller et al.l 
|2004 . All methanol maser transitions do not share the same be- 
haviour. Empirically, they form two classes (Batrlaetal. 1987). 
Class I methanol masers (e.g. at 36, 44, 84, and 95 GHz) usually 
occur in multiple locations across the star-formi ng region scattered 
around an area up to a parsec in extent (e.g ., iKurtz et al]|2004l : 
IVoronkov et alj2006l : lcvganowski et ai]|200sl) . In contrast, class II 
methanol masers (e.g. at 6.7, 12 and 107 GHz) reside in the close 
vicinity of exciting young stellar objects (YSOs) and are typically 
found as a single cl uster of emission at arcsecond resolution (e.g., 
iPhillipsetal .11 19981) . Theoretical calculations are able to explain 
this empirical classification and strongly suggest that the pumping 
process of class I masers is dominated by collisions with molecular 
hydrogen, in contra st to class II masers w hich are pumped by radia- 
tive excitation (e.g.l Voronkov et al]|2005l and references therein). 

The morphology of class I methanol masers has recently be- 
come the focus of high angular resolution studies aimed at search- 
ing for associations with other phenome na commonly observed i n 
regions of high-mass star-formation (e.g.. lCvganowski et al.l2009l) . 
The common consensus is that the majority of class I masers 
trace interface regions between outflows and molecular gas, al- 
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though direct observational evidence of thi s has been obtained 
for a limited number of sources only (e . g.. [ Plambeck & MentenI 
119901 : iKurtz et alj |2004 IVoronkov et all |2006|) . The alternative 
scenarios involving cloud-cloud collisions (e.g., ISobolevI [1993: 
ISaliietal.i r2002: Sjouwerman et al. 2010) as well as the interac- 
tion o f expanding Hll region s with the ambient molecular environ- 
ment dVoronkov et"Zll2010l) may also be realised in some sources. 
The common point of all these scenarios is the presence of shocked 
gas, where the physical conditions are favouring class I methanol 
masers (see the discussion in section |4~n below). Apart from the 
outflow associations cited above (based on the 2.12-/jm H2 emis- 
sion, which is a well known shock tracer), ICvganowski et al.l ( l2009h 
reported association of some class I maser spots with the extended 
features showing a prominent excess of the 4.5-/im emission in the 
images obtained with the Spitzer Space Telescope's Infrared Ar- 
ray Camera (IRAC), als o known as extende d green objects (EGOs ; 
ICvganowski et al.l2009h or "green fuzzies" ( [Chambers et alj2009l) . 
The excess of the 4.5-/im (IRAC band 2) emission could be a result 
of shock excitation of molecular hydrogen and carbon monoxide 
in protostellar outflows ([Cvganowski et al.l |2009[ : [Chambers et al 



[2009l:[De Buizer & Vacca|[2010l) . It is worth noting that [Chen et al 



( [2009[) demonstrated statistically the presence of an EGO in the 
vicinity of a large fraction of class I methanol masers at low an- 
gular resolution (single dish positions). To increase the number of 
class I masers studied at high angular resolution and to compare 
the morphologies observed in different maser transitions we carried 
out an interferometric survey at 36 and 44 GHz of all class I masers 
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Figure 1. Positions of class I methanol masers (red filled circles and yellow 
crosses) overlaid on top of the 3-colour Spitzer image of the G309.38— 0.13 
region. The emission in 8.0, 4.5 and 3.6 fim IRAC Spitzer bands is shown as 
red, green and blue, respectively. The filled red circles show the fitted posi- 
tion of the strongest spectral component at each location, while the weaker 
components are shown by yellow crosses (all component s are given in T a- 
ble|2). The location of a class II 6.7-GHz methanol maser ICaswelllaOOl IS 
shown by the white cross. The contours show the maximum of the 36-GHz 
emission across all spectral channels. The levels are 5, 15, 50 and 90 per 
cent of the peak 36-GHz flux density of 16.1 Jy. Low level artefacts in the 
top right corner are caused by the primary beam connection. 

previously reported in the literature and located south of declina- 
tion -35°. In this letter, we present a study of G309.38-0.13, 
an especially interesting source from the survey. It represents the 
level of morphological and kinematical complexity encountered for 
the majority of sources observed in the survey, but it has a dis- 
tinct high- veloci ty maser feature at 36 - GHz, which was quite unex- 
pected (see e.g.. lBachiller et al.ll 19901 : iKalenskii et al.ll201(]h . This 
maser is found in an area devoid of bright radio-continuum emis- 
sion, but near infrared sources which are presumably embedded 
stars. It is located approximately 20' offset from the prominent Hll 
region Gum 48d at the heliocentric distance of about 3.5 kpc (for 
further information on the region see .Karr et al.,,2009,) . 



2 OBSERVATIONS 

Observations were made with the Australia Telescope Compact Ar- 
ray (ATCA) in May 2007 as part of the interferometric survey of 
southern class I methanol masers at 36 and 44 GHz (project code 
C1642). The source was observed using the hybrid H214C array 
configuration in a few 3 minute cuts (four cuts at 36 GHz and 
six at 44 GHz) spread in hour angle for better uv-coverage. We 
did not use, for imaging, the data from the CA06 antenna (located 
about 4.5 km to the west of the other 5 antennas). The remaining 
antennas provided baselines ranging from 82 to 240 m. The posi- 
tion of the phase and pointing centre was Q2000 = 13'^47™26';01, 
feooo ~ —62° 18' 11'.' 66. We used reference pointing procedures 
and determined corrections using the continuum source 1414-59 
(which served also as a phase calibrator). From the statistics of 
pointing solutions the reference pointing accuracy was estimated 



to be 4'.'8±2'.'6. This accuracy affects the accuracy of flux density 
measurements, particularly for sources of emission which are off- 
set from the pointing centre. The positional accuracy of the maser 
locations depends on the quality of the phase calibration and is be- 
lieved to be better than 0.5 arcsec. The absolute flux density scale 
was bootstrapped from observation of Uranus (assumed flux den- 
sities were 1.38 and 1.98 Jy at 36 and 44 GHz respectively). We 
estimate the flux scale to be accurate to about 20%. The correlator 
was configured to split 8 MHz bandwidth into 1024 spectral chan- 
nels. The resulting spectral resolution and velocity range covered 
are given in Table [T] along with other observing parameters. Most 
columns of Table [T] are self-explanatory and include, in addition to 
above, the molecul ar transition and its rest frequency (adopted from 
iMuller et al.]|2004l) . full width at half maximum (FWHM) and po- 
sition angle (p. a.) of the synthesised beam reflecting the spatial res- 
olution, as well as the size of the primary beam which determines 
the field of view. The la rms noise given in the seventh column is a 
median value of rms noise levels obtained for each individual spec- 
tral plane of the image cube prior to the primary beam correction 
(i.e. with constant noise across the field of view). Scaled up with 
the primary beam, it represents the la detection limit of these ob- 
servations at a particular location. It is worth mentioning, however, 
that the ability to detect weak features is notably reduced within the 
velocity range of strong emission due to dynamic range limitations, 
which were largely a result of sparse uv-coverage. We estimate the 
dynamic range (the ratio of peak flux to rms noise in the image) to 
be about 90 and 400 for the 36 and 44-GHz image cubes, respec- 
tively. The uncertainty of the radial velocity corresponding to the 
rest frequency uncertainty is given in the ninth column of Table[T] 
Data reduction in MIRIAD followed the standard procedure in- 
cluding an opacity correction based on the model built into MIRIAD 
and atmospheric temperature, pressure and humidity registered at 
the time of observations by the ATCA weather station. We used uni- 
form weighting for the imaging. After absolute positions had been 
determined we performed self-calibration (solution interval 1 min) 
using the brightest spectral feature as a reference. The flux density 
ratio of the reference feature before and after self-calibration al- 
lowed us to estimate decorrelation caused by atmospheric phase 
variations occurring on a time-scale shorter than 5 minutes, the 
minimum temporal separation from the secondary calibrator scan. 
We estimated decorrelation factors to be l.lSiO.Ol and 1.03±0.02 
at 44 and 36 GHz, respectively. The full velocity range was imaged 
using the self-calibrated dataset. Therefore, to first order, the effect 
of the atmospheric phase stability on the derived flux densities had 
been corrected for. The residual uncertainty is expected to be small 
compared to the accuracy of the flux scale calibration. The maser 
emission was searched in the cube prior to the primary beam cor- 
rection (i.e. the noise across the field of view was constant). Then, 
the cube was divided by the primary beam model at the appropriate 
frequency and the spectra were extracted at the peak pixel by tak- 
ing a slice along the spectral axis. We followed this approach due to 
the rather high sidelobe level in the point spread function caused by 
the poor uv-coverage attained in the project. It reproduces flux den- 
sity correctly in the case of unresolved or barely resolved sources 
(i.e. smaller than the synthesised beam) and is well suited to maser 
observations. 



3 RESULTS 

The morphology of the region is shown in Figure [T] The methanol 
emission is confined to nine discrete locations, labelled A— I, 
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Table 1. Observation details. The rms noise is calculated as a median of the noise levels obtained for each spectral plane of the image cube prior to the primary 
beam correction. The velocity uncertainty corresponds to the uncertainty of the rest frequency, which is given in the brackets following the value of the rest 
frequency and is expressed in the units of the least significant figure. 
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Figure 2. Spectra of the 36 (red) and 44-GHz (blue) emission at eight locations in G309.38— 0.13 (location I is not shown here due to low signal-to-noise 
ratio in the spectrum). The dashed line shows the peak velocity of the brightest spectral feature (at location B) at the given frequency, (a) Emission from E 
(high- velocity feature), G and H was detected at 36-GHz only, (b) A, C and D were detected at 44 GHz only, (c) B (main feature) and F were detected at both 
36 and 44 GHz. The uncertainty of the relative velocity alignment between transitions is dominated by the rest frequency uncertainty of the 36-GHz transition 
shown by the horizontal error bar in spectra showing 36-GHz emission (it is comparable to the spectral resolution for the 44-GHz transition). 



spread over a 50"x30" region. The spectra at most locations are 
complex, consisting of more than one spectral feature (Fig.[2}- At 
finer scale (less than the synthesised beam size), each spectral fea- 
ture corresponds to a slightly different position in the map. Most 
of the features are believed to be spots of class I maser emission 
for reasons summarised at the end of this section. E is at essen- 
tially the same position as B, but was designated by a separate let- 
ter as it corresponds to a distinct group of velocities, blue-shifted 
by approximately 30 km s~^. We show the positions of all spec- 
tral features corresponding to each spot (circles denote the bright- 
est feature, crosses show all other features) in Fig. [T] al ong with 
the p osition of the class II methanol maser at 6.7 GHz jCaswelll 

The 6.7-GHz maser in this region peaks at —49.6 km s 
and is quite weak (I Jy). It is found at Q2000 = 13''47™23':98, 
^2000 = —62° 1 8' I2'.'0, just a few seconds of arc south of the class I 
maser emission at B and E. Both the OH and H2O masers are, 
within the measurement uncertainty, at the same position a s the 6.7- 
GHz maser and have radial velocities near —50 km s~^ jCaswelll 



ll998";'Breenetal."2010'). The 3-colour background in Fig.[T]shows 
the emission in the 8.0, 4.5 and 3.6-/im IRAC Spitzer bands as red, 
green and blue, respectively. The infrared image reveals an EGO 
in the vicinity of B and E, which is likely to be a signature of the 
shocked gas. A number of compact sources showing excess in the 
4.5-/im band (appearing green in Fig.[TJ are also present. In partic- 
ular, such a source is located near F and G. 

The spectra at all locations except I are shown in Figure[2] The 
latter is not shown due to a poor signal to noise ratio, but qualified 
as a detection due to its isolation in velocity from the other bright 
emission in the region, despite being comparable in flux density to 
dynamic range artefacts seen at other velocities. Only B and F were 
reliably detected at both 36 and 44 GHz. E, G and H were detected 
at 36 GHz only, while only 44-GHz emission was detected at A, C, 
D and I. To describe individual spectral features, we decomposed 
the spectra into a number of Gaussian components (the sum of the 
components represents each spectrum within measurement errors). 
The LSR velocity, position, velocity FWHM of the component and 
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the flux density are given on the left-hand side of Table |2] The po- 
sition fit was done assuming a point source model, typically at the 
peak velocity for the appropriate component. In the case of signifi- 
cant blending, the position fit is expected to have a systematic error 
and give a flux density-weighted average position of the blended 
Gaussian components. We estimate that the systematic errors of 
the fitted position are likely to be less than the random errors of the 
fit given in Table [2] for most components in this particular source. 
Two exceptions are the 30- Jy component towards B at 44 GHz and 
the — 78.6-km s^^ component towards E. Both were fitted with a 
single Gaussian component, although the profile shape in Figure |2] 
suggests that at least two independent components may be present 
in each case. 

The relatively compact ATCA configuration used for obser- 
vations allows us to put only a relatively weak constraint on the 
brightness temperature (> 100 K), which is insufficient to prove the 
maser nature of the observed emission unambiguously. However, 
given the line width of the majority of components (Table |2] and 
Fig. [2l(, at least some of the emission is likely to be due to maser 
action for most of the features. In addition, we verified that the 
dominant emission at both B (at both frequencies) and E is unre- 
solved at baselines to distant CA06 antenna, which implies bright- 
ness temperatures exceeding 10^ K, and, therefore, a maser origin 
of this emission. The right-hand side of Table [2] shows the param- 
eters of the line profile at each maser location, such as the LSR 
velocity and the flux density of the peak as well as the integral over 
the line profile. 



4 DISCUSSION 

4.1 High-velocity maser feature 

The blue-shifted spectral component corresponding to loca- 
tion E is a remarkable feature of this maser source. More com- 
monly, the observed velocities of class I methanol masers are 
no more than a few km s~^ from the velocity of the ambi- 
ent molecular cloud, and the velocity dispersi o n is low (e.g., 
Bachiller et al. 19 9A IPlambeck & Menteril Il990l ; IVoronkov etal] 
lOOdlArava et alJi2009ir The notable exception with a high veloc- 



ity dispersion is Gl. 6— 0.025, which has emissi on components at 
radial velocities rangin g from to 160 km s^^ jHaschick & BaanI 
ll993l : ISalii et"ai]|2002 l). However, unlike G309.38-0.13, the emis- 
sion in Gl. 6— 0.025 corresponds to distinct, but interacting, molec- 
ular clumps detect ed also in a num ber of thermal transitions of 
various molecules aSalii et al.ll2002l) . The large velocity difference 
between individual clumps reflects a violent environment and the 
exceptional kinematics in the vicinity of the Galactic centre. 

In contrast, the G309.38— 0.13 region is close to the tangen- 
tial di rection of the Sagittarius-Carina arm (see Figure 5 of lRusseill 
l2003h where the extreme velocity of E (—80 km s ^) makes an as- 
sociation with an unrelated molecular cloud highly unlikely. This 
is corroborated by the ^^CO brightness temperature limit of 0.1 K 
at the radial velocitie s below —70 km s~^ obtained with the 4- 
m Nanten telescope jSaito et alll200lh . This limit is two orders 
of magnitude lower than the peak ^^CO brightness temperature 
observed near -50 km s"^ in the G309.38-0.13 regio n (nearest 
CO clump mass is about 10"* MrT,: ISaitoetai]|200lh . which is 
close to the velocities of most of the maser spots. 

We conclude that the high-velocity 36-GHz maser in E most 
likely reflects the presence of a molecular outflow, rather than 
a separate molecular cloud projecting onto the same location by 



chance. The area around B and E is characterised by the presence 
of shocked gas traced by a strong EGO (Fig.[TJ. This supports the 
idea that an outflow is present. The close spatial location of B and E 
and the presence of the class II methanol maser (which indicates the 
presence of a high-mass YSO) in the vicinity both suggest that the 
outflow axis is close to the line of sight. 

As mentioned earlier, the interaction between outflows and the 
ambient material is believed to be responsible for the majority of 
class I methanol masers. The production of methanol by purely 
gas-phase chemical reactions is quite inefficient, but the passage of 
weak shoc ks is expected to rele ase methanol from the grain man- 
tles (e.g.. ^Hart quist et^ll995h . Enhanced methanol abundance 
arising from shocks driven by outflows into the ambient cloud 
has also been con firmed observationally (e.g.. lGibb & Davisil993 : 
lGaravetaI.|[200"a) . In addition, outflows compress the medium, 
which makes collisions more frequent and the pumping mech- 
anism more efficient. However, methanol is expected to survive 
sputtering and desorption of grain mantles only if the shocks 
are relative ly mild, i.e. at s hock velocities not greatly exceeding 
10 km s"^ jGarav et al.ll200a) . Therefore, the detection of a maser 
feature offset by 30 km s~^ with respect to the velocity of the am- 
bient gas is somewhat surprising. 

This discrepancy can be resolved if the outflow interacts with 
a moving parcel of gas. In this case, the same radial velocity of 
the maser can be achieved for a lower shock velocity. The mecha- 
nism of this gas acceleration is unknown at present. It could be by 
a previous outflow event or a different outflow from another YSO. 
In addition, the clustered environment of high-mass star formation 
opens a range of possibilities for gravitational acceleration of gas 
clumps (e.g. mechanisms which caused radial velocity outliers ob - 
served in clusters of low-mass protostars by ICovevetal.ll200^ . 
One may speculate that the curved morphology of the EGO as- 
sociated with B and E (Fig.[TJ suggests interaction with an outflow 
coming from the south-west (the general direction of the Hll region 
Gum 48d). This may be the additional outflow responsible for the 
unusual kinem atics. However , the distance of Gum 48d is estimated 
to be 3.5 kpc jXarr et alj|2009h . and therefore Gum 48d must be at 
least 20 pc away from the maser. A separation of the terminal shock 
from the driving s ource as large as thi s is unprecedented (e.g. the 
outflow studied bv lBrooks et al.l ( 120031 ) extends only up to 1.5 pc) 
and suggests that, if the EGO is indeed shaped by the second out- 
flow, the driving source is not the star responsible for the excitation 
of Gum 48d. 



4.2 Relative flux density of 36 and 44-GHz methanol masers 

An inspection of Table |2] and Fig.[2]suggests that the 36 to 44-GHz 
flux density ratio varies significantly from component to compo- 
nent. The majority of maser features were detected in a single tran- 
sition only and most striking is the lack of a 36-GHz counterpart to 
the brightest 44-GHz spectral feature (Fig|2};). Although the veloc- 
ities corresponding to E (high-velocity feature) were not covered 
by the 44-GHz measurement, the preliminary results of a follow 
up observation carried out with the Mopra 22-m single dish in Jan- 
uary 2010 showed no 44-GHz methanol (as well as no 48-GHz 
methanol, SiO 1-0, and CS 1-0) emission near —80 km s~^ above 

1.3 Jy. With the caveat about possible temporal variability over the 
2.5-year timescale, this situation represents another extreme case 
of a 36-G Hz maser without a 44-GHz counterpart. The numerical 
models o f lSobolevetai] ( l2005l) suggest that the 44-GHz masers are 
quite sensitive to the geometry of the maser region. Unlike 36-GHz 
masers, which can be produced in the spherically symmetric case. 
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Table 2. Fit results and profile parameters. The uncertainties are given in parentheses and expressed in units of the least significant figure. Notes: (a) the 
uncertainty is half of that for the line FWHM, (6) the uncertainty is the spectral resolution listed in Table[T] 
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-49.62 


20.7 


(2) 


li 


!06 


(1) 


1.1 


(2) 


0.57 


(8) 
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4-1 - 3o E 


36 


-60.56 


19.27 


(5) 


li 


! 11.5 


(3) 


1.98 


(8) 


1.54 


(6) 


-55.16 


1.9(4) 


7 


(1) 








-55.29 


19.32 


(4) 


li 


! 11.6 


(3) 


2.2 


(1) 


1.31 


(5) 
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36 


-51.50 


21.24 


(3) 


17 55.0 


(2) 


2.04 


(5) 


1.48 


(3) 


-50.39 


3.3 (2) 


3.79 


(5) 








-51.432 


21.52 


(2) 


17 56.9 


(1) 


0.275 (8) 


3.30 


(9) 










1 


7o - 6i A+ 


44 


-53.25 


25.19 


(5) 


li 


! 10.8 


(2) 


0.81 


(8) 


0.43 


(4) 


-53.29 


0.6 (2) 


0.404 (8) 



an optical depth along the line of sight higher than in the orthogonal 
direction seems to be the key property required to produce bright 
44-GHz masers. One could expect this elongation of the gas clumps 
involved in the maser action to be somehow related to the outflow 
orientation and morphology. However, the complexity of the maser 
spot distribution (see Fig.[TJ does not allow us to interpret the flux 
density ratios for each spot in a consistent manner based on just 
the spot location with respect to the suspected outflow (the position 
of the 6.7-GHz maser presumably marks the location of the YSO 
driving the outflow roughly along the line of sight). 



5 CONCLUSIONS 

(i) Spots of class I methanol maser emission in G309.38— 0. 13 
are spread over an area 50" X 30" in extent (see Tableland Fig.©. 
Both 36 and 44-GHz maser emission was detected at two locations, 
while the emission from the remaining locations was detected at a 
single frequency only (see Fig.|2j. 

(ii) We report the detection of a high-velocity spectral feature 
at 36 GHz, blue-shifted by about 30 km s^^ from the peak veloc- 
ity. This is the largest velocity offset reported so far for a class I 
methanol maser source associated with a single molecular cloud. 
The maser corresponding to this feature (labelled E) is located 
within a few seconds of arc of the maser at B (the brightest spectral 
feature at both frequencies) indicating a possible association with 
an outflow parallel to the line of sight. 

(iii) The results highlight the importance of not only high spec- 
tral and spatial resolution for studies of class I methanol masers, 
but also ample velocity coverage and a field of view wider than 



the typical primary beam size of available instruments. All obser- 
vational data gathered so far lack at least one (and sometimes all) 
of these properties. In addition, the results show that the 36 and 44- 
GHz transitions are very complementary to each other. The high- 
velocity spectral feature in G309.38— 0. 13 is a good example of a 
phenomenon which has been missed by previous maser surveys. 
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